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ABSTRACT 

The Universe's Epoch of Reionization (EoR) can be studied using a number of observational 
probes that provide complementary or corroborating information. Each of these probes suffers 
from its own systematic and statistical uncertainties. It is therefore useful to consider the mu- 
tual information that these data sets contain. In this paper we present a cross-correlation study 
between the kinetic Sunyaev-Zel'dovich effect (kSZ) - produced by the scattering of CMB 
photons off free electrons produced during the reionization process - and the cosmological 
21 cm signal - which reflects the neutral hydrogen content of the Universe, as a function 
of redshift. The study is carried out using a simulated reionization history in 100 IrT 1 Mpc 
scale N-body simulations with radiative transfer. In essence we find that the two probes anti- 
correlate. The significance of the anti-correlation signal depends on the extent of the reioniza- 
tion process, wherein extended histories result in a much stronger signal compared to instanta- 
neous cases. Unfortunately however, once the primary CMB fluctuations are included into our 
simulation they serve as a source of large correlated noise that renders the cross-correlation 
signal insignificant, regardless of the reionization scenario. 

Key words: cosmology: theory, cosmic microwave background, large scale structure of Uni- 
verse, diffuse radiation, radio lines: general 



1 INTRODUCTION 

The Epoch of Reionization (EoR) is one of the least explored pe- 
riods in the history of the Universe. At present, there are only a 
few tentative observational constrains on the EoR like the Gunn- 
Peterson troughs ( |Gunn & Peterson| [19651 |Fan et al.|[2006| ) and 
the cosmic microwave background (CMB) E-mode polarization 
(Page et al. 2007) at large scales. Both of these observations pro- 
vide strong yet limited constraints on the EoR. In the near future, 
however, a number of observations at various wavelengths (e.g., 
redshifted 21 cm from HI, Lyman-a emitters, high redshift QSOs, 
etc.) are expected to probe this pivotal epoch in much greater de- 
tail. Among these, the cosmological 21 cm transition line of neutral 
hydrogen is the most promising probe of the intergalactic medium 
(IGM) during reionization ( Ma dau et al.|1997| > 

A number of radio telescopes (e.g. LOFAbQ MWAQ and 
SKA0 are currently being constructed/designed that aim at detect- 
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ing the redshifted 21 cm line to study the EoR. Unfortunately, these 
experiments will suffer from a high degree of contamination, due 
both to astrophysical interlopers like the Galactic and extra-galactic 
foregrounds and non-astrophysical instrumental effects (e.g. |Jelic| 
|et al.||2008[ |Labropoulos et al.||2009| >. Fortunately, the signal has 
some characteristics which differentiates it from the foregrounds 
and noise, and using proper statistics may extract signatures of 
reionization (e.g. Furlane tto et al.|2 004, Harker et al. 2009b a). In 
order to reliably detect the cosmological signal from the observed 
data, it is essential to understand in detail all aspects of the data and 
their influence on the extracted signal. 

Given the challenges and uncertainities involved in measur- 
ing the redshifted 21-cm signal from the EoR, it is vital to cor- 
roborate this result with other probes of the EoR. In this paper 
we study the information imprinted on the CMB by the EoR and 
its cross correlation with the 21 cm probe. Given the launch of 
the recent PLANCK satellite, that will measure the CMB with un- 
precented accuracy, it is only fitting to conduct a rigorous study into 
the cross-correlation of these data sets. 

One of the leading sources of secondary anisotropy in the 
CMB is due to the scattering of CMB photons off free electrons, 
created during the reionization process, along the line of sight 
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(LOS) (Zeldovich & Sunyaev 1969 ). These anisotropics when in- 
duced by thermal motions of free electrons is called the thermal 
Sunyaev-Zel'dovich effect (tSZ) and whence due to bulk motion 
of free electrons, the kinetic Sunyaev-Zel'dovich effect (kSZ). The 
latter is far more dominant during reionization (for a review of sec- 
ondary CMB effects see |Aghanim et al.|2 008). 

The kSZ effect from a homogeneously ionized medium, i.e., 
with ionized fraction only a function of redshift, has been stud- 
ied both analytically and numerically by a number of authors; the 
linear regime of this effect first calculated by Suny aev & Zel 



dovich (1970 ) and subsequently revisited by Ostriker & Vishniac 



|dov 

( 1986J and |Vish niac (1987) - hence also referred to as the Ostriker- 
Vishniac (OV) effect. In recent years various groups have calcu- 
lated this effect in its non-linear regime using semi-analytical mod- 
els and numerical simulations ( Gnedin & Jaffe 2001 ; Santos et al. 
2003 ; Zhang et al. 2004). These studies show the contribution due 
to non-linear effects being important only at small angular scales 
(/ > 1000), while the OV effect dominate at large angular scales. 

The kSZ effect from patchy reionization was first esti- 
mated using simplified semi- analytical models ( [Santos et al.|2 003 ) 
wherein they concluded that fluctuations caused by patchy reion- 
ization dominate over anisotropics induced by homogeneous reion- 
ization. However, for a complete picture of the CMB anisotropics 
induced by the EoR a more detailed modeling is required. Over 
and above the underlying density and velocity fields these details 
should include the formation history and "nature" of the first ioniz- 
ing sources and the radiative transport of ionizing photons to derive 
the reionization history (sizes and distribution of the ionized bub- 
bles). Some recent numerical simulations of the kSZ effect during 
the EoR were carried out by |Salvaterra et al.] ( |2005| >; |Zahn et aT] 
^j^ ' ^o^^^i^ ' ^ssv et al.H2007^ 

Cross-correlation between the cosmological 21 cm signal and 
the secondary CMB anisotropics provide a potentially useful statis- 
tic. The cross-correlation has the advantage that the measured 
statistic is less sensitive to contaminants such as the foregrounds, 
systematics and noise, in comparison to "auto-correlation" studies. 
Analytical cross-correlation studies between the CMB temperature 
anisotropics and the EoR signal on large scales (I ~ 100) have been 
carried out by |Alvarez et al.| ( |2006j>;|Adshead & Furlanettolf2 008) 
and on small scales (/ > 1000) by[Cooray (2004 ); Salvaterra et al.| 
( 2005); Slosar et al. (2007). Thus far the only numerical study of the 
cross-correlation has been carried out by Salvaterra et al. (2005). 
Some additional analytical work on cross-correlation between the 
E- and B -modes of CMB polarization with the redshifted 21 cm sig- 
nal have been done by Tashiro et al. (2008 ); Dvorkin et al. (2009 ). 

In this paper we first calculate the kSZ anisotropics from the 
homogeneous and patchy reionization based on 100 /i -1 Mpc scale 
numerical simulations of reionization. We then cross-correlate 
them with the expected EoR maps obtained via the same simula- 
tions, and we discuss how the large-scale velocites and primary 
CMB fluctuations influence the cross-correlation. Although simi- 
lar in some aspects, the work presented here differ from Salvaterra 
|et al.| ( [2005 1 ) substantially. First, Salvaterra et al. used a relatively 
small computational box (20 Mpc) incapable of capturing rel- 
evant large-scale density and velocity pertubations. Secondly, the 
primary CMB fluctuation, that manifests itself as a large back- 
ground noise, is not taken into the cross-correlation study. And fi- 
nally, there is a difference in procedure for calculating that cross- 
correlation coefficient. 

The paper is organized as follows. In Section [2] we discuss 
kSZ signal and cosmological 21 cm signal from the EoR. In Sec- 
tion [3] we present the numerical simulations employed to obtain 



the kSZ and EoR maps for a specifc reionization history. Cross- 
correlation between the cosmological 21 cm fluctuations (EoR sig- 
nal) and the kSZ anisotropics, together with the influence of the 
large-scale velocities and the primary CMB fluctuations on CMB- 
EoR cross-correlation are discussed in Section [4] Finally in Sec- 
tion [5] we present our discussions and conclusions on the topic. 

Thro ughout we assume A CDM-cosmology with WMAP5 pa- 
rameters ( [Komatsu etal.|2 008 ): H = 71.9 kms _1 Mpc _1 , Q b = 
0.0441, Q m = 0.258 and Q A = 0.742. 



2 THEORY 

In the following section we briefly review the theoretical aspects 
of the kinetic Sunayev-Zel'dovich (kSZ) effect and the cosmolog- 
ical 21 cm signal from the epoch of reionization. In addition, the 
relevant mathematical forms used to calculate the kSZ and the cos- 
mological 21 cm signals are presented. 



2.1 Kinetic Sunayev-Zel'dovich effect 

The temperature fluctuation of the CMB caused by the Thompson 
scattering of its photons off populations of free electrons in bulk 
motion, for a given line of sight (LOS), is: 



ST 
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r n e (r ■ v)dt, 
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where r is the optical depth of electrons to Thomson scattering, v 
the bulk velocity of free electrons and r the unit vector denoting 
the direction of the LOS. The integral is performed for each LOS 
with t T the time at the epoch recombination and to the age of the 
Universe today. Note that all quantities are in physical units. Tem- 
perature fluctuations produced at time t will be attenuated due to 
multiple scattering along the LOS to the present time and is ac- 
counted for by the e _T term. 

The electron density can be written as the product of the neu- 
tral atom density n n and ionization fraction x e . Both n n and x e 
vary about their average value n n and x e and thus these fluctua- 
tions can be defined as S = n n /n n — 1 and S Xe = x e /x e — 1 and 
consequently the electron density expressed as: 



e(l + 6 + 6 Xe +55 x 



(2) 



In first approximation one can just follow the reionization of 
hydrogen and assume that the neutral atom density equals neutral 
hydrogen density. However, in our simulation we follow both hy- 
drogen and helium reionization. Assuming both hydrogen and he- 
lium follow the underlying dark matter density, the neutral atom 
density is a sum of the neutral hydrogen (nm) and neutral helium 
(riHei) densities: n n = (mn+miei)(l + £). Moreover, the electron 
density can be written as: 



TlmXmi + ^Hel^Hell + 2n H eI^HeIII, 



(3) 



where 7iHii,Heii,Hein and XHii,Heii,Hein are densities and ioniza- 
tion fractions of HII, Hell and Helll respectively. The ionization 
fractions are defined as: xhii = ^hii/^hi, ^Heii = ^Heii/^Hei 

andXHelll = TlHelll/nHel. 

The mean hydrogen and helium densities vary with redshift 

as n H i,Hei = n H i(o),Hei(o)(l + z) 3 , where n H i(o),Hei(o) are the 
mean hydrogen and helium densities at the present time: n H i(o) = 
1.9 x 10~ 7 cm -3 andnHei(o) = 1.5 x 10~ 8 cm -3 . 
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Figure 1. A slice through the simulated redshift cube of the the cosmological 21 cm signal (top panel) and the kSZ effect (bottom panel) in the case of the 
'Stars' patchy reionization model. The angular scale of the slices is ~ 0.6°. 



By inserting Eq. [2] into Eq. [T] and converting Eq. [T] from an 
integral in time to one in redshift space ^\(z), we get: 



ST 
T 



p (1 + z) 2 _ T _ 
= -cr T n n(0 ) j — — — e x e 

•(1 + 5 + 5 Xe + 65 Xe )v r dz, 



(4) 



3 + Qa where H is the present 
is the matter and the dark 



where v T is component of v along the LOS (v r = r • v) and n n ( ) = 
^hi(o) + ^Hei(o)- For a ACDM Universe the Hubble constant at 
redshift z is H = H \/Q m (l + z 
value of the Hubble constant, Q u 
energy densities, respectively. 

For homogeneous reionization histories, i.e. a uniform change 
in the ionization fraction as a function of redshift, Eq.|4]becomes: 

ST\ _ f z ° (1 + z^ 2 

— - I -cr T n H (o) / 



kSZ 



-CTTn H (o) 



- r x e (l + ^) Vr dz, (5) 



which means that the kSZ fluctuations are induced only by spatial 
variations of the density field. The linear regime of this effect is 
called the Ostriker-Vishniac (OV) effect. The OV effect is of second 
order and peaks at small angular scales (arc minutes) and has an 
rms of the order of a few /jK. 



dt -- 



In order to make transformation of the Eq.[T]to the redshift space we use 



dz 

~H{z)[l+z] 



, where H (z) is Hubble constant at redshift z. 



2.2 The Cosmological 21 cm signal 

In radio astronomy, where the Rayleigh- Jeans law is applicable, 
the radiation intensity, I(y) is expressed in terms of the brightness 
temperature X^: 



/(!/) = 



kT b , 



(6) 



where v is the frequency, c is the speed of light and k is Boltz- 
mann's constant. The predicted differential brightness temperature 
of the cosmological 21 cm signal with the CMB as the background 
is given by ( |Field| 1958| [T959] |Ciardi & Madau|2003| ): 



ST h = 
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Here T s is the spin temperature, xm is the neutral hydro- 
gen fraction, S is the matter density contrast and h — 
Ho I (100 kms _1 Mpc _1 ). If we express the neutral hydrogen frac- 
tion as xhi = xhi(1 + <Wi)> Eq.|7]becomes: 



6 Tb = 26mKx m (l + S + S x 
Tomb A / Qbh 2 



1 
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In his two seminal papers, [Field| (T958 1959 ) calculated the 
spin temperature, T s , as a weighted average of the CMB, kinetic 
and colour temperatures: 



TcMB + VkinTkin + VaTcx 
1 + Vkin + Vol 



(9) 



where Tcmb is the CMB temperature and ytin and y a are the 
kinetic and Lyman-a coupling terms, respectively. We have as- 
sumed that the color temperature, T a , is equal to Tki n jMadauet al.| 
1997). The kinetic coupling term increases with the kinetic temper- 
ature, whereas the y a coupling term depends on Lyman-a pumping 
through the so-called Wouthuysen-Field effect (Wouthuysen 1952, 
Field 1958). The two coupling terms are dominant under different 
conditions and in principle could be used to distinguish between 
ionization sources, e.g., between first stars, for which Lyman-a 
pumping is dominant, vs. first mini-quasars for which X-ray pho- 
tons and therefore heating is dominant (see e.g., |Madau et al.|1997[ 
|Zaroubi et al.|2007||Thomas & Zaroubi|2008| ). 



3 SIMULATIONS 

The kSZ (5T/T) and the cosmological 21 cm maps (5Tb) 
are simulated using the following data cubes: density (5), ra- 
dial velocity (v r ) and HI, HII, Hel, Hell and Helll fractions 
(#Hi,Hii,Hei,Heii&Heni). The data cubes are produced using the 
BEARS algorithm, a fast algorithm to simulate the EoR signal 
( [Thomas et al.|2009) . 

In the following subsections we summarize the BEARS algo- 
rithm and describe operations done on the output in order to cal- 
culate the kSZ and the EoR maps. Furthermore, we detail the cal- 
culations to obtain the optical depth and kSZ signal along a certain 
LOS. Finally we present the maps of the kSZ temperature fluctu- 
ations for the two patchy reionization models ('stars' and 'mini- 
quasars') and discuss aspects of their contribution to the signal. 



3.1 BEARS algorithm: overview 

BEARS is a fast algorithm to simulate the underlying cosmologi- 
cal 21 cm signal from the EoR. It is implemented by using an N- 
body/SPH simulation in conjunction with a 1-D radiative transfer 
code under the assumption of spherical symmetry of the ionized 
bubbles. The basic steps of the algorithm are as follows: first, a cat- 
alogue of ID ionization profiles of all atomic hydrogen and helium 
species and the temperature profile that surround the source is cal- 
culated for different types of ionizing sources with varying masses, 
luminosities at different redshifts. Subsequently, photon rates em- 
anating from dark matter haloes, identified in the N-body simula- 
tion, are calculated semi-analytically. Finally, given the spectrum, 
luminosity and the density around the source, a spherical ioniza- 
tion bubble is embedded around the source, whose radial profile is 
selected from the catalogue as generated above. For more details 
refer |Thomas et al.| ( [2009l >. 

As outputs we obtain data cubes (2D slices along the fre- 
quency/redshift direction) of density (5), radial velocity (v r ) and 
hydrogen and helium fractions (xHi,Hii,Hei,Heii&Heiii)- Each data 
cube consists of about 850 slices each representing a certain red- 
shift between 6 and 11.5. Slices have a size of 100 h -1 comoving 
Mpc and are defined on a 512 2 grid. Becauses these slices are pro- 
duced to simulate a mock dataset for radio-interferometric experi- 
ments, they are uniformly spaced in frequency (therefore, not uni- 
form in redshift). Thus, the frequency resolution of the instrument 



dictates the scales over which structures in the Universe are aver- 
aged/smoothed along the redshift direction. The relation between 
frequency v and redshift space z is given by: 

V2\ 



1, 



(10) 



where vn — 1420 MHz is the rest frequency that corresponds to 
the 21cm line. 

The final data cubes are produced using approximately 35 
snapshots of the cosmological simulations. Since choice of the red- 
shift direction in each box is arbitrary, three final data cubes can be 
produced in this manner (x, y and z). 



3.2 Randomization of the structures 

The kSZ effect is an integrated effect and is sensitive to the struc- 
ture distribution along the LOS. To avoid unnatural amplification 
of the kSZ fluctuations due to repeating structures in the simulated 
data cubes, we follow |Iliev et aL] ( |2007| ) and introduce random- 
ization of the structures along the LOS over 100 Mpc/h scale in 
two steps. First, each 100 Mpc/h chunk of the data cube is ran- 
domly shifted (assuming periodic boundary conditions) and rotated 
in a direction perpendicular to the LOS. The shift can be positive 
or negative in any direction (x and(or) y) by an integer value be- 
tween and 512. While the rotation can be clockwise or anticlock- 
wise by an rnv/2 angle (n=0, 1,2,3). Secondly the final data cube is 
produced by assembling the first 100 Mpc/h part from the x-data 
cube, second from the y-data cube, third from the z-data cube and 
then back to the x-data cube and so on., to a distance that spans the 
comoving radial distance between redshifts 6 and 12. 



3.3 Optical depth 

Thomson optical depth r at redshift z is: 



C(7T 
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JO 



H(z) 



-dz, 
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where c = 2.998 x 10 8 m s _1 is the speed of light, ctt = 
6.65 x 10 -29 m 2 the Thomson scattering cross section for elec- 
trons, n e the density of free electrons and H(z) the Hubble con- 
stant at redshift z. 

In our simulations we split the integral into two parts. The first 
part that represents the mean Thomson optical depth (foe) between 
the redshift 0-6 and the second, tq z , from redshift 6 to a desired 
redshift z. This choice is driven by the limited redshift range (z ~ 
6 — 11.5) of imminent radio astronomical projects designed to map 
the EoR. Under the assumption that the reionization is completed 
by the redshift 6, the mean Thomson optical depth f 06 is 0.0517. 
Note that our simulation is set to have a mean Thomson optical 
depth of 0.087, as obta ined from the CMB data (r = 0.087=b0.017, 
|Komatsu et"aL] ( [2008l >). 



3.4 Creating the kSZ and EoR maps 

For clarity we summarize the steps we follow to create the kSZ and 
EoR maps for a given scenario of the reionization history: 

(i) Using the output of BEARS, data cubes for the density, radial 
velocity, helium and hydrogen fractions are produced. 

(ii) Data cubes are randomized over 100 Mpc/h scale along the 
redshift direction. 
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Figure 2. The same as a Fig.^but for the 'QSOs' patchy reionization model. 
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(iii) Using Eq.[TT]the Thomson optical depth, r, is calculated to 
a redshift z. 

(iv) Using the integrand of the Eq. [4] data cubes with the kSZ 
signal is produced as a function of redshift. 

(v) Integrating along each LOS through the kSZ data cube the 
integrated kSZ map is obtained. Note that we assume that the reion- 
ization is complete by redshift 6, so the integral in Eq. [4] spans 
z > 6. 

(vi) Finally, the brightness temperature fluctuations, 5Tb, is cal- 
culated using Eq.[8] 

In the following sections we will use the kSZ and EoR maps 
produced for three different models of reionization: 

(i) HOMOGENEOUS: Reionization history is homogeneous and 
ionized fraction follows: 

Xe = k(l-z ■ ( lT} 

1 ~\~ e { 2 reionJ 

with Zreion set to 8.5 and k = 1, 2, ... tunes the "rapidness" of the 
reionization process. 

(ii) PATCHY Stars: Reionization history is patchy, gradual and 
extended with stars as the sources of ionization. 

(iii) PATCHY QSOs: Reionization history is patchy and rela- 
tively fast with QSOs as the ionizing sources. 

Fig.[T]& Fig.|2]show slices through the simulated redshift cube 
of the cosmological 21 cm signal (STh) and the kSZ effect (£Tksz) 
in the case of 'Stars' and 'QSOs' patchy reionization models. The 
angular size of the slices is ~ 0.6°. 



Apart from the difference in the global shape of the reioniza- 
tion histories driven by 'Stars' and 'QSOs' (see Fig. [4]), the average 
size of the ionization bubble is also smaller in 'Stars' compared 
to that of 'QSOs'. For a detailed description and comparison of 
reionization histories due to 'Stars' and 'QSOs' see Thomas et al. 
( [20091 ). 

The kSZ anisotropics from patchy reionization is induced by 
both fluctuations of the density field 5 and ionization fraction S Xe 
(see Eq. [4|. [Santos et alT] < [2003| > found that kSZ anisotropics from 
5 Xe fluctuations dominate over 5 modulated fluctuations (OV ef- 
fect). In order to test this result with our simulations we split the 
integral in Eq.[4]into three parts and produce three integrated kSZ 
maps (for the 'Stars' model see Fig.[3j. The first term '1 + 5' repre- 
sents the density induced secondary anisotropics (OV effect). The 
'Sxe' term represents the secondary anisotropics due to patchiness 
in the reionization and '££ Xe ' represents a higher order anisotropy. 

The mean and rms of the '1 + 5\ 'Sx^ and '£(5^' com- 
ponents of the simulated kSZ maps are shown in the Table [T] for 
patchy reionization in the 'Stars' and 'QSOs' model. The rms 
value of the maps is used as a measure of the fluctuations. We con- 
firm that the '(5^' fluctuations are indeed larger than density in- 
duced anisotropies ('(5') for both patchy reionization models. How- 
ever the difference between the '^ e ' and '£' fluctuations is much 
larger for the 'Stars 'reionization history model than 'QSOs' model. 
Also note that the third order anisotropy CSSx^) is not negligible 
in both reionization scenarios. For completeness we also show con- 
tribution from the pure Doppler term ('!') in Eq.[4] 
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Figure 3. The simulated kSZ anisotropics induced by '1 + 5' (first panel), '5 Xe ' (second panel) and 'Sd Xe ' term (third panel) in Eq.^for patchy reionization 
model by 'Stars'. The kSZ anisotropics induced by all terms together in Eq.^are showed on fourth panel ( 'TOTAL '). The mean and rms of the simulated kSZ 
maps is shown in TableU\ Note that each map has its own color scale. 



Table 1. The mean and rms of the '1+5', i 5 Xe ' and i 55 Xe ' simulated kSZ 
maps for both 'Stars' (see Fig. [3j and 'QSOs' patchy reionization model. 
Co is a cross-correlation coefficent at a zero lag between coresponding kSZ 
maps and integrated EoR map (see Section |4j. For completeness we also 
show results for the pure Doppler term ('1') in Eq.[4] 
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Figure 4. The mean ionization fraction x e as a function of redshift for the 
'Stars' and 'QSOs' patchy reionization model. 



4 CROSS-CORRELATION KSZ-EOR MAPS 

The kSZ effect from the EoR is expected to be correlated with cos- 
mological 21 cm maps for a homogeneous reionization history and 
anti-correlated wh en patchy (|Coora y 2004, Salvat erra et al.| 2005 ; 
|Alvarez et all2006l|Slosar et al.|2007||Adshead & Furlanetto|2 008 ). 
In this section simulations described in Section |3] are used to ex- 
plore the small angular scale cross-correlation between the kSZ 
effect and EoR maps for five different reionization histories. Fur- 
ther, we will fold-in the influence of i) the large-scale velocities on 
the kSZ effect and ii) the primary CMB fluctuations on the cross- 
correlation. 

Throughout the paper we will use a normalized cross- 
correlation in order to be able to compare results from a different 
pairs of maps. The normalized cross-correlation between two im- 
ages (aij and bij) with the same total number of pixels n is defined 
at zero lag as: 

c = ^-Y Kj-")^-^ (13) 

Tl - 1 ^— ' CFaCfb 

where a (b) is the mean and a a (<Jb) the standard deviation of the 
image a (b). However, the cross-correlation between the kSZ and 
the EoR map needs to be considered more carefully. 

The fluctuations of the kSZ effect over the simulated map are 
both positive and negative, since the radial velocity v r can be both 
positive and negative (see Eq. |4j. In contrast, the EoR signal fluc- 
tuations in our simulations are always positive (see Eq. [8]). When 
calculating the cross-correlation between these two maps, we are 
interested in finding the number of points at which both signals are 
present (homogeneous reionization model) or where one signal is 
present and the other absent (patchy reionization model). In other 
words only the absolute value of the kSZ fluctuation is relevant in 
our calculation and not its sign. 

4.1 Homogeneous reionization history 

We explore the cross-correlation between the kSZ map and inte- 
grated EoR map in the case of three different homogeneous reion- 
ization histories (HRH1, HRH2 & HRH3). These histories are 
modulated by Eq. [12] with k controlling the duration of reion- 
ization. The mean ionization fraction x e (z) for these models are 
shown in Fig. [5] 

The cross correlation between an integrated kSZ map and an 
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Figure 5. The ionization fraction x e as a function of redshift for three dif- 
ferent models of the homogeneous reionization history (HRH 1, HRH 2 & 
HRH 3). All three models are defined by Eq.^H^but have different values of 
k (different reionization durations). The half of the reionization process is 
set in all three models to z re i on =8.5. 

integrated EoR map results in a coefficient C HRH j = 0.10 ± 
0.03 for an extended homogeneous reionization history (HRH 1). 
For HRH 2 C 0jHRH 2 = 0.21 ± 0.02 and HRH 3 C Q HRH 3 = 
0.24 ± 0.02. The errors are estimated by performing a Monte Carlo 
calculation with 200 independent realizations of the integrated kSZ 
and EoR maps using the randomization procedure explained in Sec- 
tion|3l 

As expected, the integrated kSZ and EoR maps are correlated 
for homogeneous models of reionization. Furthermore, the corre- 
lation depends on the duration of reionization with larger values 
for more 'rapid' reionization. These results are in agreement with 
|Alvarez et aE] ( [2006l >. 

4.2 Patchy reionization history 

For patchy reionization models we first cross-correlate the kSZ and 
the EoR map at a given redshift. The resulting zero lag coeffi- 
cient (Co), as a function of redshift, is shown in Fig. [6] The solid 
black line represents the correlation for 'Stars' while the dashed 
red line the 'QSOs' patchy reionization model. As expected for 
patchy reionization in both models, the kSZ and the EoR map anti- 
correlate at individual redshifts. 

Note that Eq.[l3]is not taking into account a statistical weight 
(significance) of the result. The statistical weight can be defined 
with a factor 2n*/n where n* represents the number of aij • by- 
products that are non-zero. This makes the result the most signifi- 
cant when both signals are represented equally and less significant 
when one signal is present and the other not. 

The obtained anti-correlation is also evident by visual inspec- 
tion of the kSZ and EoR slices through the simulated redshift cubes 
(see Fig.[T]&[2]). One can see that the kSZ signal is present only at 
the regions where the EoR signal is not. This result is not surprising 
since the EoR signal is proportional to neutral hydrogen while the 
kSZ to the ionized, both of which are almost mutually exclusive. 



Figure 6. The zero lag cross-correlation coefficient (Co) between the kSZ 
map and the EoR map at a given redshift. The solid black line corresponds 
to the 'Stars' while the dashed red line the 'QSOs' patch reionization model. 
For both reionization models we find an anti- correlation between the maps. 

In reality, we are not able to measure the kSZ effect at a cer- 
tain redshift but only an integrated effect along the entire history. 
Thus we can only cross-correlate the integrated kSZ map with the 
integrated EoR map and/or the EoR maps at different redshifts^] 

Fig. [7] shows the integrated EoR and kSZ map for the 'Stars' 
(first two panels) and 'QSOs' (last two panels) patchy reionization 
model . The cross-correlation coefficient at zero lag for these two 
models are Co, stars = —0.17 and Co,qso s = —0.02. In order 
to determine the error on the kSZ-EoR cross-correlation, we per- 
form a Monte Carlo calculation. After creating 200 independent 
realizations of the integrated kSZ and EoR maps using the random- 
ization procedure explained in Section [3] we calculate the cross- 
correlation coefficient for each pair of realizations. Finally, we cal- 
culate the mean and standard deviation of the cross-correlations. 
For the 'Stars' model we get Co, stars = — 0.16 =b 0.02, while for 
the 'QSOs' model C ,qso s = -0.05 ± 0.02. 

To understand higher values of the cross-correlation coeffi- 
cient in 'Stars' compared to the 'QSOs' model, one needs ana- 
lyze Fig. [4] and Table [T] It is evident from Fig. [4] that the reioniza- 
tion history is gradual and extended with stars as ionizing sources, 
compared to a shorter and sharper history with QSOs as ioniz- 
ing sources. Moreover, the patchy term C5 Xe ') of the kSZ fluc- 
tuations is much larger than the homogeneous component in the 
anisotropy ('£') in the case of 'Stars' than for 'QSOs' model (see 
Table [TJ. We showed beforehand that the kSZ effect correlates with 
the cosmological 21 cm signal for homogeneous reionization and 
that the correlation is strongest for an 'instant' reionization history. 
We also obtain the same result by correlating different kSZ com- 
ponents with the integrated EoR map (see Table [TJ. Combining 
these results we see the cross-correlation is driven by the patchy 
kSZ anisotropics in the 'Stars' model, while in the 'QSOs' model 
the homogeneous and patchy kSZ anisotropics tend to cancel each 

5 This is because, unlike the kSZ effect, we can potentially obtain redshift- 
specific information of neutral hydrogen via upcoming radio telescopes. 
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Figure 7. The integrated EoR and kSZ map for the 'Stars' (first two panels) and 'QSOs' patchy reionization model (second two panels). The mean cross- 
correlation coefficient at the zero lag between integrated EoR map and integrated kSZ map is Co, stars — —0.16 zb 0.02 for the 'Stars' and Co,qsOs — 
-0.05 ± 0.02/or the 'QSOs' model 



other. As a consequence, the anti-correlation in 'QSOs' model is 
much weaker than that of 'Stars'. 

In addition to the balance between homogeneous and patchy 
kSZ anisotropics that governs the (anti-)correlation between the 
kSZ and the EoR maps, the size of the ionized bubbles also play 
a key role. Recall that the average size of the ionization bubble 
is larger for 'QSOs'. As a result, the underlying structure within 
the ionized bubble will additionally reduce the anti-correlation and 
might change the scale of (anti-)correlation. 

From now on we will just concentrate on cross-correlations 
using 'Stars' since the 'QSOs' model does not show a significant 
anti-correlation. Fig. [8] shows the correlation coefficient as a func- 
tion of lag (C(0)) between the integrated kSZ and the integrated 
EoR map. The dashed red lines represents the estimated error ob- 
tained from Monte Carlo simulations. As in Salvaterra et al. (2005), 
we find that the two signals are anti-correlated below a characteris- 
tic angular scale C and this scale indicates the average size of the 
ionized bubble which in our case is 6 C ~ 10 arcmin. 

Salvaterra et al. ( 2005) also showed that the amplitude of the 
anti-correlation signal increases with decreasing redshift and that 
the characteristic angular scale shows a redshift evolution. In order 
to test this in our simulation, we calculate the redshift evolution 
of the zero lag cross-correlation coefficient between the integrated 
kSZ map and the EoR map at different redshifts (Fig. [9). To cal- 
culate the error in the cross-correlation, we generate 200 differ- 
ent realizations of the kSZ and corresponding EoR cubes using the 
randomization procedure explained in Section [3] Then, around a 
desired redshift we fix the kSZ effect to zero and integrate along 
non-zero part of the kSZ cube. Finally, we cross-correlate the in- 
tegrated kSZ map with the EoR map at the desired redshift and 
estimate the error on the cross-correlation between the integrated 
kSZ map and the EoR map at the certain redshift. 

From Fig. [9] we find no coherent redshift evolution of the anti- 
correlation signal and on the contrary at few redshifts the two sig- 
nals correlate instead of anti-correlating. The correlation at a given 
redshift is caused by i) the patchy nature of the EoR signal, which 
implies that there are some redshifts at which the EoR map contains 
none or only a few small ionized bubbles. If one correlates such an 
EoR map with the integrated kSZ map, the outcome is a correlation 
between the two, and because of an insignificant number of the 
ionization bubbles there is no contribution to the anti-correlation, 
ii) the kSZ signal being patchy. There are some redshifts where 
the kSZ signal from a certain ionization bubble does not contribute 
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Figure 8. The cross-correlation between integrated EoR and integrated kSZ 
map as a function of lag ( C (0) )for the 'Stars ' reionization history scenario 
(dashed white line). The gray shaded surface represents the estimated error 
obtained by Monte Carlo simulation. Note that the correlation coefficient at 
the zero lag is C = -0.16 ± 0.02. 



significantly or at all to the integrated kSZ map (see Fig.[l0|. This 
could happen due to a weak kSZ signal from a certain ionized bub- 
ble or due to cancellation of the kSZ signal from another ionization 
bubble along the LOS. 

The analysis is repeated for different binnings in frequency of 
the EoR map but the results does not differ significantly. We also 
calculate the redshift evolution of the charateristic angular scale 
(# c ), but we do not find any coherent evolution. This result is driven 
by the fact that the contribution of the kSZ signal from a certain 
redshift to the integrated kSZ map is not significant or is even non- 
existent. As a result, if there is no coherent redshift contribution to 
the integrated kSZ map there will be no coherent redshift evolution 
of the kSZ-EoR cross-correlation signal (see Fig.[T0|l. 

The discrepancy between our results and that of Salvaterra 
|et al.| p005| > is mainly due to the difference in the method to cal- 
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Figure 9. The redshift evolution of the zero leg correlation coefficient be- 
tween the integrated kSZ map and the EoR map at the certain redshift. The 
result is shown for the 'Stars' reionization history model. Note that the EoR 
map at the certain redshift is produced by integrating 100 h~ 1 Mpc vol- 
ume around the desired redshift. 



culate the cross-correlation coefficient. In Salvaterra et al., they 
first calculate the cross-correlation coefficient (not normalized) be- 
tween a certain kSZ and EoR map. Then, they scramble both maps 
without keeping any structural information and calculate the cross- 
correlation coefficient. They compare the coefficients in the two 
cases to draw their conclusion. In contrast to Salvaterra et al., 
we first calculate the normalized cross-correlation coefficient (see 
Eq. between a pair of kSZ-EoR map. And then for compar- 
ison, we preform a Monte Carlo simulation to generate different 
realizations of the kSZ and the EoR maps . However, despite the 
cross-correlation procedure used, once the primary CMB fluctua- 
tions are included we are not able to find any significant kSZ-EoR 
cross-correlation (see Sec.|4.4|). 



4.3 Large-scale velocity 

Our simulation volume is (100 Mpc) 3 (see Section^. Thus, 
large-scale velocities associated with bulk motions, on scales > 
100 h' 1 Mpc are missing. The missing velocities represent ~ 50% 
of the total power in the velocity field as given by the linear theory. 

|Iliev et aL] ( 2007| ) showed that the large-scale velocities on 
scales > 100 h~ Mpc increases the kSZ signal. Motivated by this 
result, we approximately account for the missing large-scale veloc- 
ities as follows: first, we assume that every 100 /i -1 Mpc chunk 
of our simulation cube has a random large-scale velocity compo- 
nent vls- Since our simulation cube is produced using 15 simula- 
tion boxes (100 /i -1 Mpc), we need in total 15 vls- We randomly 
choose a realization of the 15 ^ls based on a velocity field power 
spectra from linear theory. By doing this we ensure that the veloc- 
ities are correlated at large scales. Finally we add the missing vls 
component to each 100 h~ x Mpc chunk of the simulated cube. 

Based on 200 realizations of the large-scale velocity field, we 
have found that the large-scale velocities increase the kSZ signal 
during the EoR by 10%. But on average we do not find any signif- 
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Figure 11. The map of the primary CMB fluctuations generated as a Gaus- 
sian random filed with a power spectrum obtained from the CMBFAST al- 
gorithm. 



icant increase or decrease in the kSZ-EoR cross-correlation. How- 
ever, for ~20% of all large-scale velocity realizations we find an 
increase in the cross-correlation signal by a factor two or larger and 
for ^2% a factor three or larger. 



4.4 Primary CMB 

Up to now, in our cross-correlation analysis we only considered 
secondary CMB anisotropics generated by the kSZ effect. In the 
actual experiment, the CMB data will comprise not only the kSZ 
anisotropies which are secondary, but also the primary and other 
secondary CMB anisotropies (for a recent review see Aghanim 
|et al.| < [2008] >). In this subsection we will examine the influence of 
the primary CMB fluctuations on the detectability of the kSZ-EoR 
cross-correlation. 

We simulate the primary CMB fluctuations in the following 
way: first the CMB power spectra is obtained using CMBFAST 
(U. Seljak & M. Zaldarriaga, 2003) and then the map of the pri- 
mary anisotropy is produced as a random Gaussian field with this 
power spectrum. An example of the simulated primary CMB map 
is shown on the Fig.[TT] The size of the map coresponds to the size 
of the simulated EoR and kSZ maps. Notice the leak in power at 
small scales due to Silk damping (Sil k|1967) . 

In order to calculate the noise in the cross-correlation intro- 
duced by the primary CMB fluctuations, we generate 200 dif- 
ferent realizations of the primary CMB (pCMB) fluctuations. We 
then add secondary kSZ anisotropies induced by the 'Stars' (map 
shown in the Fig. [TJ, and calculate the cross-correlation between 
the pCMB+kSZ map and the integrated cosmological 21 cm map. 
The obtained zero lag cross-correlation coefficient is 0.0 =b 0.3. 
The noise introduced by the primary CMB fluctuations is too large 
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Figure 10. TOP PANEL: A random line of sight through 'Stars' kSZ cube, ST^qz( z )> that is averaged over 10 pixels (~ 0.7 arcminj at each redshift. BOTTOM 
PANEL: For the same LOS the cumulative integral of the kSZ effect, J* #?kSZ (z)dz . Note that there is no coherent redshift contribution to the integrated 
kSZ effect. 

to find any significant kSZ-EoR (anti-)correlation. However, one 
has to remember that the primary CMB anisotropics are damped 
on small angular scales and that on these scales the secondary 
anisotropics are the dominant component of the CMB power spec- 
tra (see Fig. [12]). Utilizing this fact, one can do a cross power spec- 
trum and see the correlation as a function of angular scale. Pursuing 
this lead, we calculate the kSZ-EoR cross spectrum first without 
and then with the primary CMB added to the kSZ map. 

The cross spectrum (C*) between the two images of a small 
angular size is give by: 

Cf ~ P fc X = — A ™ ' ( 14 ) 

p,q£k 

where A p , q is Fourier transform of the first image, B* q the com- 
plex conjugate of the Fourier transform of the second image and 
rtk is number of points in the k-th bin (k = ^/p 2 + q 2 ). Note 
that we assume a 'flat-sky' approximation (e.g. [White et al. 1999 ): 
k 2 P(k) ^J^rCi \i=2irk which is valid for I > 60. 

Fig. [13] shows the cross power spectrum between the kSZ 
anisotropics and the integrated cosmological 21 cm map for reion- 
ization due to 'Stars'. It is evident from the plot that the two im- 
ages anti-correlate on large scales (/ < 8000) but that the anti- 
correlation becomes weaker towards the smaller angular scales. At 
angular scales I > 8000, there is no (anti-)correlation. 

We also calculate the cross power spectrum between inte- 
grated EoR map and integrated kSZ map with primary CMB fluc- 
tuations included. However, the noise introduced by the primary 
CMB is too large to find any significant correlation at scales I < 
8000. 

This result might be driven by the simulation box size and 
reionization scenarios considered in this study and does not mean 
that a cross-correlation signal is absent at all scales and reionization 
histories. In order to test this, one needs to explore the kSZ-EoR 
cross-correlation using simulation with the box size that is larger 
than the one used in this study (100 /i -1 Mpc). 
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Figure 12. The power spectra of primary CMB fluctuations (dotted line) 
and kSZ anisotropics obtained from the simulated maps (solid line). 

In the following subsection we try to use different cross- 
correlation methods and filtering techniques to explore the kSZ- 
EoR cross-correlation at scales where kSZ signal dominates over 
primary CMB fuctuations and where correlation signal is non-zero 
(4000 < I < 8000). 



4.5 Additional cross-correlation techniques 

For better understanding of the properties of the kSZ-EoR cross- 
correlation, and with the hope of being able to find the cross- 
correlation signal in the presence of primary CMB fluctuation, in 
this subsection we apply techniques of filtering, wavelet decompo- 
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Figure 13. The cross spectrum (see Eq. \l4\ between integrated kSZ map 
and integrated EoR map for the 'Stars' reionization history (dashed white 
line). The gray shaded surface represents the estimated error obtained by 
Monte Carlo simulation. Note that the primary CMB fluctuations are not 
included. 



sition and relative entropy on our data. We will only use the inte- 
grated kSZ map and integrated EoR map from the 'Stars' model of 
reionization, since the strongest cross-correlation signal is obtained 
for this model in the previous sections. Note that in a following 
analysis we first use only the kSZ and the EoR map and then as a 
second step we include the primary CMB fluctuations. 

Fig. [14] shows the zero lag cross-correlation coefficient as a 
function of the three different filtering procedures. The first one 
uses a high-pass, the second a low-pass and the third a band-pass 
filter that passes out only a certain scale. In all three cases the filter 
is based on the 'Top hat' function. We filter out the desired scale 
from both the kSZ map and the EoR map and calculate the cross- 
correlation coefficient at zero lag. The results are shown for the 
low-pass and high-pass filter as a function of the FWHM of the 
filter and for the band-pass filter as a function of scale. 

The plot on the left panel in Fig. [14] implies that the anti- 
correlation is strongest on the largest scales of the map. By adding 
smaller scales, the correlation coefficient decreases meaning that 
smaller scales introduce noise in the correlation. The middle panel 
in Fig. [14] suggests the same behaviour. By removing the large 
scales, the cross-correlation signal becomes quite small. Finally, 
the third panel of Fig. [14] suggests that the large scales are indeed 
the dominant component of the anti-correlation signal. 

As a next step in our analysis we include the primary CMB 
fluctuations. However we obtain the same result as discussed in 
the previous subsection. On the scales where the kSZ anisotropics 
dominate over primary, or the anti-correlation signal is too weak 
or the noise introduced by residuals of the primary CMB fluctua- 
tions is too large to find any statistically significant kSZ-EoR (anti- 
correlation. 

The wavelet analysis of the maps is done using Daubechies 
and Coiflet wavelet functions. Both integrated kSZ map with added 
primary CMB flucutations and integrated EoR map are decom- 
posed to a certain wavelet mode and then they are cross-correlated. 



Because the outcome is similar to that of filtering we will not dis- 
cuss this further. 

The last method applied to the data is the "relative entropy", 
also know as Kullback-Leibler distance. The relative entropy is a 
measure of the shared information between two variables (two im- 
ages) by comparing the normalized distribution of the two. This 
method too did not produce any significant result. 



5 DISCUSSION AND CONCLUSIONS 

This paper presents a cross-correlation study between the kinetic 
Sunyaev-Zel'dovich (kSZ) effect and cosmological 21 cm signal 
produced during the epoch of reionization (EoR). The study uses 
an N-body/SPH simulation along with a 1-D radiative transfer code 
(the BEARS algorithm, [Thomas et al.| ( |2009| >) to simulate the EoR 
and to obtain maps of the cosmological 21 cm signal and of the 
kSZ effect. The maps are produced using the 100 h' 1 Mpc co- 
moving simulation box for five different (3 homogeneous and 2 
patchy) models of reionization history. The homogeneous model 
with varying degree of "rapidness" of the reionization process is 
given by Eq. [12] The patchy reionization histories include one by 
'Stars' (gradual) and the other by 'QSOs' (instant). 

For a homogeneous reionization history we find that the kSZ 
map and the integrated EoR map are correlated. Furthermore, that 
the correlation depends on duration of reionization with larger val- 
ues for more "rapid" models. This result agrees with the analytical 
kSZ-EoR cross-correlation analysis carried out by | Alvarez et al.| 
( [20061 ). 

For patchy reionization models we find that the kSZ temper- 
ature fluctuations are of the few \iK level (see Table [TJ, and is 
in agreement with previously obtained results by Salvaterra et al. 
(2005|; [lliev et al.| ( [2007] ). In addition, we show that the tempera- 
ture fluctuations induced by the patchiness of the reionization pro- 
cess C$x e ' term in Eq.[4]) is larger than the density induced fluctua- 
tions (homogeneous '1 + 5' term in Eq.[4]). The difference between 
the two is stronger for the extended history ('Stars' model) than in 
more rapid reionization histories ('QSOs' model) (see Table[T]). 

As a first step in the kSZ-EoR cross-correlation study of 
patchy reionization histories we cross-correlate the kSZ map and 
EoR map at each redshift (see Fig. [T] &[2j. As expected, the kSZ 
and the EoR map anti-correlate at certain redshifts (see Fig. [6). 

We then cross-correlated the integrated cosmological 21 cm 
map and the integrated kSZ map for patchy reionization (see 
Fig. [7]). The two signals show significant anti-correlation only in 
the 'Stars' model (Co, stars = -0.16 ± 0.02, C ,qso s = -0.05 ± 
0.02.). The result is driven by the balance between homogeneous 
and patchy ('l+£' and '£ Xe ' terminEq.[4]) kSZ anisotropies and the 
average size of the ionized bubbles. Since the homogeneous kSZ 
anisotropies correlate and patchy kSZ anisotropies anti-correlate 
with the cosmological 21 cm maps, the two effects tend to can- 
cel each other. In addition the average size of the ionization bub- 
ble is larger for 'QSOs' than in 'Stars' model and the structure of 
matter within the ionized bubble reduce the cross-correlation. As 
a consequence the kSZ-EoR anti-correlation is much stronger for 
the extended ('Stars' model) reionization history than for a more 
instant history ('QSOs' model). 

For a patchy model of reionization we estimated the redshift 
evolution of the correlation coefficient (Co) and characteristic an- 
gular scale 0c. This was done by cross-correlating the integrated 
kSZ maps with the EoR maps at different redshifts (see Fig. [9}. In 
contrast to Salvaterra et al. (2005), we do not find any significant 
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Figure 14. The zero lag cross-correlation coefficient as a function of the three different filtering procedures. The first one uses a high-pass filter, the second one 
uses low-pass filter and the third a band-pass filter that passes only a certain scale. In all three cases the filter is based on the 'Top haf function. The dashed 
white line is mean and the gray shaded surface represents the estimated error obtained by Monte Carlo simulation. 



coherent redshift evolution of Co and 6c . The discrepancy between 
the results is caused by the difference in a procedure used for cal- 
culating cross-correlation. However, despite the cross-correlation 
procedure used, once the primary CMB fluctuations are included 
we are not able to find any significant kSZ-EoR cross-correlation. 

The influence of the missing large-scale velocities on the kSZ 
signal and kSZ-EoR cross-correlation was investigated. Although, 
the large-scale velocities increase the kSZ signal by 10%, we do 
not find, on average, any significant change in the kSZ-EoR cross- 
correlation. However, for ~20% of large-scale velocity realizations 
we find an increase in the cross-correlation signal by a factor two 
or larger and for ^2% a factor three or larger. 

Data from CMB experiments contains both the secondary (e.g. 
kSZ) and primary anisotropies. For completeness of our study we 
calculated the noise in the kSZ-EoR cross-correlation introduced 
by the primary CMB fluctuations and found that its addition re- 
duces the cross-correlation signal to zero (Co = 0.0 ± 0.3). The 
cross-correlation was also performed on scales where the kSZ 
anisotropies dominate over the primary CMB fluctuations (/ > 
4000, see Fi g.[T2] l. We have done this by calculating cross-power 
spectra (Fig. |13|>, applying different filtering methods (Fig. [T4j on 
the data and by doing wavelet decomposition. However, the out- 
come of the analysis is that on the scales where the kSZ anisotropies 
dominate over primary, or the anti-correlation signal is too weak 
or the noise introduced by residuals of the primary CMB fluctua- 
tions is too large to find any statistically significant kSZ-EoR (anti- 
correlation. 

As a sanity check we calculate the kSZ-EoR cross-correlation 
using the simulation obtained by |Iliev et ah] ([2007 ) ('f250C 
100 /i -1 Mpc simulation). The reionization history of this model 
is similar to our 'QSOs' model. The reionization history is rela- 
tively sharp and instant. The cross-correlation coefficient at zero 
lag for the integrated kSZ map and integrated EoR map is Co = 
—0.04 ± 0.02. The result is in agreement with the result obtained 
from the 'QSOs' model. We also calculated the redshift evolution 
of the zero lag cross-correlation coefficient and have found no co- 
herent redshift evolution. 

In view of all the results obtained from the kSZ-EoR cross- 
correlation study, we conclude that the kSZ-EoR anti-correlation 
on scales captured by our simulation box (~ 0.6°) is not a reliable 
technique for probing the EoR. However, there is still hope that we 



will be able to find the correlation between the kSZ and EoR signal 
on scales larger than ~ 1°, where the patchiness of the ionization 
bubbles should average out (see Alvarez et al. (2006 ) and Tashiro 
et al., in preparation). Finally, it is important to note that the kSZ 
signal induced during the EoR could still be detected in the power 
spectra of the CMB and used to place some additional constrains 
on this epoch in the history of our Universe. 
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